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A  STUDY  OF  SIMPLE  WAVES 
J.  Hardy 

Lawrence  Radiation  Laboratory^  University  of  California 
Livermore,  California 

August  10,  1965 

ABSTRACT 

A  method  for  construction  of  simple  waves  is  presented.  Attention  is  focused  on 
forward -facing  compression  waves,  specifically  at  the  points  at  which  shocks  develop, 
and  what  to  expect  concerning  shock  strength,  etc.  Discussion  is  limited  to  ^-law  gases, 
although  the  technique  is  general. 
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INTRODUCTION 

Section  1  gives  the  construction  of  simple  wave  solutions  for  y -law  gases.  The 
method  is  general;  restriction  to  y-law  gases  is  for  convenience.  Velocity,  density,  .  .  .  , 
are  constant  along  a  family  characteristics  presented.  The  region  multiply  covered  by  the 
family  of  characteristics  is  discussed  along  with  the  envelope  of  the  family. 

Section  2  presents  the  solution  (complete  with  graphs)  to  several  piston  motions, 
illustrating  the  properties  of  the  most  general  solutions. 

Section  3  shows  that  shocks  develop  only  on  the  envelope  of  the  characteristics. 

The  local  behavior  of  the  envelope  is  given,  depending  on  the  sign  of  a  certain  differential 
inequality.  The  behavior  is  examined  specifically  for  a  neighborhood  of  zero  time. 

Means  are  presented  for  analysis  of  a  general  flow. 


SECTION  1 

The  hydrodynamic  equations  for  one -dimensional  flow  are 
Pt  +  up^  +  pu^  =  0. 

where  p(x,  t)  is  the  density,  u(x,  t)  is  the  material  speed,  e(x,  t)  is  the  specific  energy 
density,  and  p(x,  t)  is  the  pressure. 

In  the  case  of  isentropic  flow,  the  system  yields  p  =  p(p);  and  for  a  y-law  gas, 

p  =  •  Under  the  enumerated  circumstances 

^\,PqJ 


where  c 


2 


p^  +  up 
+UU 
=  7p/p. 


X 

X 


+  pu^  =  0. 


(2) 
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The  interest  is  in  solutions  to  this  latter  set  of  equations.  The  flow  considered 
will  be  that  of  a  piston  pushing  into  a  gas  at  rest,  the  piston  speed  g(t)  being  such  that 
the  flow  is  everywhere  continuous. 

The  general  character  of  these  solutions  is  well  known.  They  are  simple  waves, 
flows  whose  image  in  the  hodograph  plane  (p,  u)  is  a  curve.  If  g(t)  is  such  that  the  piston 
starts  from  zero  speed  (necessary)  and  moves  to  the  right,  the  following  relations  among 
u,  c,  p,  p,  Uq,  Cq,  Pq,  Pq  are  known  to  hold  (u^,  p^,  c^,  p^  are  the  rest  values): 


P 


U  -  Uq  \  2y/ (7-1) 


P 


’  ^0 


1  + 


-  1 


u  -  Uq\2/(7-1) 

J 


(3) 


relations  which  hold  for  a  foreward -facing  simple  wave. 

Since  the  value  of  Uq  makes  no  difference  in  any  arguments,  we  shall  assume 
Uq  =  0.  Substituting  the  expression  for  p  as  a  function  of  u  into  the  continuity  equation 

Pt  +  up^  +  pu^  =  0, 


we  obtain 


o' 


a  quasi-linear  partial  differential  equation  whose  solution  u(x,  t)  is  constant  along  the 
characteristic  curves 


1 


dx 

TJJ. 

2 


U 


These  are  straight  lines,  since  u  is  constant  along  them. 

Consider  now  the  initial  value  problem  (see  Fig.  1).  Att  =  0,  u  =  0  along  the  x  axis. 
Along  x(t),  u  =  g(t),  the  piston  speed.  Along  t  =  0,  the  characteristics  are  the  lines 
X  =  Xq  +  c^t,  and  along  these  lines  u,  the  solution  to  equation  (4)  is  zero. 


At  time  t^,  the  piston  is  at  x(tQ)  =  j  g(?  )  dS  with  speed  g(tQ).  The  character¬ 
istic  passing  through  this  point  is  ^ 

X  =  x(tQ)  +  I^Cq  +  (t  -  Iq),  (5) 

and  along  it  u  =  g(tQ).  Thus,  a  solution  has  been  constructed,  a  ruled  surface,  giving 
u(x,  t)  for  those  points  in  the  (x,  t)  plane  for  which  the  simple  wave  is  a  valid  solution. 
Since  near  t  =  0,  g(t)  >  0,  the  one-parameter  family  of  characteristics  constructed  does 
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Fig.  1.  Construction  of  characteristics. 

not  cover  the  plane  simply.  Some  points  in  the  (x,  t)  plane  must  be  covered  at  least 
twice.  It  is  the  implication  of  this  we  wish  to  investigate.  Note  that  if  g  is  discontinuous, 
double-valued  points  (x,  t)  exist  beyond  this  time.  The  double-valued  region  for  u(x,  t) 
gives  physically  unacceptable  solutions.  At  some  points  of  the  boundary  of  this  region, 
a  shock  will  develop.  The  simple  wave  ceases  to  be  a  valid  solution  to  the  system  (1) 
in  the  whole  (x,  t)  plane.  ^  ' 

Let  the  piston  be  located  at  Xj(t)  at  the  timet,  x^(t)  =  ^  g(5  )  d| ,  g(?  )  is  monotone 
increasing.  The  characteristic  passing  through  (t,  x^(t))  intersects  the  characteristic 
through  (0,  0),  the  wave  front  traveling  with  velocity  c^,  at  a  time 

'  (r+)"grt)  [^^^0 

Let  t  tend  to  zero.  If  g’  (0)  f  0,  then 


"  iy  +  l)g'  (0)  ■ 

Otherwise,  the  demand  that  g{t)  >0,  t  >  0,  requires  that 
lim  T,  (t)  =  +  00  . 


The  characteristics  passing  through  (t,  x^(t))  and  (f,  Xj^(t'  ))  i 


intersect  at  a  time 


T^a,  t') 


2(t  -  t') 


The  envelope  (r^  ,  X2)  of  such  lines  is  given  by  the  limit  as  t‘  approaches  t, 

^2<*>  'TfTTJgTt)  h  *  8<‘>  * •«'  <*>] 


tg'  (t) 
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where  g'(t)  f  0,  or,  if  g'  (t)  =  0, 

lim  T  n(t)  =  +<)o.  (10) 

t'-^t  ^ 

The  coordinate  x  for  the  envelope  is  obtained  from  substituting  T2(t)  in  equation  (5). 

For  the  class  of  monotone  increasing  functions  g(t)  (giving  rise  to  forward-facing 
compression  waves)  the  time  at  which  a  shock  appears  in  the  flow,  a  time  corresponding 

to  inf(t,  x)  for  all  (t,  x)  in  the  region  doubly  covered  by  the  characteristics,  is  given  by 
(t) 

the  time 

T„  =  inf  T„(t).  (11) 


x„  =  infx(T„)  located  by  equation  (5).  For  purposes  of  examining  these  functions,  the 

<^2> 

following  easily  developed  formulae  are  useful: 


d  /.V  --y  -  1  .  0 

dt  “7  +  1  (7  +  l)  g(t) 


.m. 

g(t)^ 


- 


where  tj^  t  in  equation  (14). 


SECTION  2 

Prior  to  considering  general  classes  of  problems,  we  will  consider  several 
examples . 

Example  1:  c^  =  1,  g(t)  =  2t,  7  =  2. 

The  solution  to  the  problem  is  given  by 

Jo,  X 

u(x,  t)  = 

[2tg,  X  =  tQ*^  +  (1  +  3tQ)(t  -  y,  tQ>  0, 

p  and  p  being  given  by  equation  (3),  and  it  solves  the  piston  problem  in  the  region  of 
single  valuedness.  This  solution  is  plotted  in  Fig.  2.  The  piston  path,  some  character¬ 
istics,  and  the  point  where  the  shock  originates  are  indicated.  Suppose  t'  >  t  >  0. 
Formula  (8)  gives  the  time  of  intersection  for  the  characteristics  passing  through  the 
piston  at  these  times  as 
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u 


Fig.  2.  Example  1. 


T^a.v )  =  |<t  + 1')  +  ^  >  I . 

with  equality  occurring  for  t  =  t*  =0.  Thus,  any  characteristic  passing  through  a  point 
near  the  origin  intersects  the  wave  front  at  a  time  r  1/3  and  at  a  point  x(t)  1/3, 
these  values  being  the  limiting  values.  For  x>l/3,  t>l/3,  u(x,  t)  is  double -valued  on 
the  wave  front.  At  t  =  1/3,  a  shock  starts  at  x  =  1/  3,  with  jump  in  velocity  growing 
continuously  from  zero. 

Example  1  is  a  special  case  of  g(t)  =  2at,  a  >  0,  where  always  the  infimum  of  the 
times  of  intersection  of  pairs  of  characteristics  is  r  -  the  shock  always  originating 
on  the  wave  front.  Note  that  this  point  recedes  to  infinity  as  a  tends  to  zero. 

Example  2:  c^  =  1,  7  =  3,  g(t)  =  (1  -  2t)  -  1.  The  characteristics  all  pass 

through  the  point  (l/2,  1/2).  Also,  lim  g(t)  =  +00  .  This  solution  is  graphed  in  Fig. 

t^l/2-0 

3.  The  example  ceases  to  be  of  interest  at  t  =  1/2. 

From  it,  an  interesting  example  can  be  constructed.  Let 

■  3  |(1  -  -  1,  0  <  t<  1/4, 

•-^y2  - 1,  t  >  1/4. 
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Fig.  3,  Example  2. 


For  t  >1/4,  the  piston  follows  a  straight  line  in  the  (x,  t)  plane,  and  the  characteristics 
are  all  parallel  straight  lines,  intersecting  x  =  c^t  at  times  at  least  1/4,  Along  these 
u(x,  t)  has  the  value  ^f2  ~  1,  The  solution  thus  consists  of  this  plane  parallel  to  the  (x,  t) 
plane,  plus  that  section  of  example  2  corresponding  to  characteristics  for  time  t  <  1/4. 

At  t  =  1/2,  a  pair  of  shocks  appear  at  x  =  1/2,  one  with  jump  in  velocity  smaller 
than  ^J2  -  1  but  still  positive.  This  shock  proceeds  into  the  undisturbed  medium  with  this 
initial  shock  strength.  Characteristically,  a  weaker  shock  proceeds  from  this  point  back 
into  the  disturbed  medium  -  toward  the  piston.  A  contact  discontinuity  is  present  always 
at  the  point  of  origination  of  these  two  shocks.  Until  the  weak  shock  reaches  the  piston, 
the  shock  velocities  are  constant,  the  states  behind  each  being  constant. 

Example  2  is  an  example  of  a  piston  velocity  with  an  important  separation  property 
to  be  discussed  directly. 
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Example  3:  =  1,  7  =  3,  and  g(t)  defined  by 


g(t)  = 


|^4t  -  1.05, 


0  ^  t  <  0,3, 
t  >  0.3, 


a  continuous  function  with  jump  in  (t)  at  t  =  0,3.  The  solution  is  graphed  in  Fig.  4 
complete  with  characteristics  and  piston  path.  Examination  of  the  characteristics 
through  the  formulae  (8),  (12),  and  (13)  shows  that  the  double- valued  region  in  the  range 
of  u  has  a  point  on  its  boundary  at  t  =  0.44375  .  .  .  and  x  =  0.20937  .  .  ,  ,  all  other  points 
on  its  boundary  (and  interior)  corresponding  to  larger  times.  This  point  is  on  the 
characteristic  passing  through  the  piston  location  at  time  t  =  0.3;  u(x,  t)  has  a  discontinu¬ 
ous  derivative  at  the  piston  at  t  =  0.3  and  this  discontinuity  propagates  along  this  charac¬ 
teristic. 


A  shock  originates  at  this  boundary  point  with  shock  strength  (jump  in  velocity, 
etc.  )  growing  continuously  from  zero.  Observe  that  this  shock  originates  in  the  region 
of  continuous  flow  well  behind  the  wave  front. 

From  subsequent  discussion  it  will  be  apparent  that  it  would  be  possible  to  modify 
g(t),  t  >  0.3,  in  such  a  fashion  that  a  shock  (see  Example  2)  appears  at  this  point 
t  =  0.44375,  x  =  0.20937,  with  a  nonzero  jump  in  velocity, 

Another  important  observation  is  that  the  shock  can  originate  in  the  region  of  con¬ 
tinuous  nonconstant  flow  (behind  the  wave  front)  without  a  discontinuity  in  any  derivative 
of  g(t)  occurring.  In  example  3,  take 
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g(t)  =  t  >  0. 

The  boundary  point  of  the  region  of  double  values  with  minimum  time  occurs  at  t  =  0.97847 
.  .  .  ,  X  =  0.9642  .  .  .  ,  which  is  a  point  in  the  region  of  continuous  nonconstant  flow. 

These  examples  exhibit  the  important  properties  of  the  solutions  to  piston  problems 
which  persist  for  a  time  as  simple  waves,  with  shocks  then  forming  somewhere.  Be¬ 
cause  of  the  arbitrariness  of  g(t)  and  the  possibility  of  patching  functions  of  dilferent 
characters  together, flows  with  quite  varied  behavior  are  possible,  as  for  example, 
simultaneous  appearance  at  t  =  tg  of  n  shocks  of  varying  initial  shock  strengths,  etc. 

(sum  of  the  jumps  of  velocity  at  the  shocks  is  strictly  less  than  g(t)). 

Before  discussing  local  behavior  of  the  envelope,  let  us  show  that  all  candidates 
for  infima  are  provided  by  the  envelope  of  characteristics. 


SECTION  3 


Consider  Fig.  5,  and  piecewise  smooth  g(t),  g’(t>  f  0.  Let  C  and  C  be  a  pair  of 
characteristics  corresponding  to  times  t  and  t'  .  The  slope  of  the  characteristics  is  a 
monotone  increasing  function  of  this  time,  the  time  they  pass  through  the  piston  location, 


Thus,  for  the  time  t^  =  (t  +  t*  )/  2  the  slope 
is  intermediate.  It  must  intersect  both  C 
and  C  at  times  that  may  be  equal  but  can- be 
no  later  than  the  time  corresponding  to  the 
intersection  of  C  and  C’  .  Choose  the 
earlier  time.  Pick  a  time  equidistant 
from  t^  and  the  time  corresponding  to  the 
characteristic  giving  the  earliest  time. 

The  characteristic  C^  corresponding  to  t^ 
then  intersects  C^  and  either  C  or  C  (which¬ 
ever  gave  the  least  time)  at  a  time  no  later 
than  the  time  of  intersection  of  this  later 


Fig.  5.  pair.  Thus,  continuing,  we  construct  a 

convergent  sequence  t^,  ^  ^n '  *  *  ' 

the  sequence  of  whose  pairs  of  characteristics  is  converging,  and  the  sequence  of  whose 
points  of  intersection  is  converging  to  a  point  on  the  envelope  of  the  characteristics. 
Thus,  for  any  such  pair  C,  C'  of  characteristics,  no  matter  how  remote,  there  is  a 
point  on  the  envelope  of  the  characteristics  whose  time  is  not  larger  than  the  time  of 
intersection  of  the  pair. 

The  time  T^(t)  is  the  time  of  intersection  of  the  characteristic  corresponding  to 
t  with  (Fig.  6);  hence,  for  each  time,  by  the  above,  there  is  a  value  of  r  (t)  which  is 
no  larger.  Thus,  the  infimum  sought  over  the  whole  field  of  intersection  characteristics 
will  be  furnished  by  inf  r^(t). 

4.  " 


-9- 


Fig.  6. 

Turn  now  to  a  discussion  of  local  behavior.  Suppose  all  characteristics  corres¬ 
ponding  to  a  time  Iq  >  0  pass  through  the  same  point  (r,  x(t)).  The  function  g(t)  from 
equation  (9)  must  satisfy  the  differential  equation 

■:^+^^g(t)  +  (t-T)g'(t)  =  o 

which  has  a  solution  for  tQ  <  t  <  r  , 


This  solution  satisfies  the  equation 


o' (  I  r  -  t  J 


jXzl 

7+1 


7-1 


_  iL 


g 


g(t) 


(16) 


(17) 
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As  it  happens,  much  can  be  learned  concerning  the  local  behavior  of  intersecting 
characteristics  (a  neighborhood  of  characteristics  about  that  corresponding  to  t^)  accord¬ 
ing  as 

^  [‘^0 

is  positive,  negative,  or  zero.  If  it  is  positive,  in  every  sufficiently  small  neighbor¬ 
hood  of  Iq  to  the  left  of  t^.  g(t)  <  g^(t).  If  it  is  negative,  the  inequality  is  reversed. 

This  property  enables  a  discussion  of  the  local  behavior  of  the  envelope  of  characteristics 
to  be  given. 

Consider  Fig.  6.  Suppose  that  in  a  neighborhood  of  l^,  g(t)  is  twice  continuously 
differentiable  and  the  expression  (18)  is  positive.  Then,  from  formulae  (13),  (14)  and 
(15)  it  is  clear  that  all  of 


dT2(t) 

~dt 


dx(T2) 

dt 


,  and 


are  positive.  Thus,  as  t  increases  from  t^,  x(t2)  and  increase  monotonically. 
Since  ^  <  x  (72)  “  ^  '^2 '  increasing,  the  point  (T2,  x  (t^))  is  moving  to 

the  left  of  Cg,  that  is,  into  the  doited  region.  Thus  72(1^)  furnishes  a  minimum  for 
72(1^),  those  values  of  t  in  the  neighborJiood  of  t^^  which  are  larger  than  tj^. 

If  the  expression  (18)  is  negative,  by  the  same  reasoning,  the  point  is  moving  into 


the  crosshatched  region  (to  the  right  of  C2  with  both  7p  and  x(72)  decreasing  monotoni¬ 
cally  as  t  increases  from  tj^). 


In  either  case,  the 
of  these  cliaractei'islics 


curve  72.  tangent  to  C2  at  t^,  since  it  is  the  envelope 

it  IS  also  easily  verified  from  equations  (13)  and  (14)  that 


dx(7„) 

-A — —  =  + 

d72  0 


Observe  that  a  local  maximum  or  minimum  in  ^2(1)  corresponds  to  a  cusp  in  the 
envelope. 

We  will  now  consider  a  neighborhood  of  the  origin  in  which  g(t)  is  twice  continuously 
differentiable  (permitting  g"(t)-*oo  as  t-  0)  and  g'  can  be  zero  only  at  the  origin. 

Observe  that  if  g'(0)  ^  0.  7^(0)  =  72(0)  -  2Cq/ (y  +  l)g'  (0).  If  g'  (0)  ^  0.  -+00  as 

t->  0  but  remains  finite  for  all  subsequent  values  of  g(t).  72(1)  tends  to  infinity  as  t 
tends  to  any  value  of  for  which  g'(t)  is  zero. 

If  g'  (0)  >  0,  then,  according  as  tiie  quantity  (18)  is  positive  or  negative,  72(1), 
x(72)  moves  to  the  left  or  right  of  C^,  starting  from  the  point  2c  J (y  +  1)  g'  (0).  2Cq 
/{y  h  l)g'  (0).  In  the  first  case  botii  7,,(t)  and  7j^(t)  increase  and  a  minimun  is  furnished 
by  T2(0).  In  the  second  case,  72(1)  is  rnoiiotone  decreasing,  so  the  infimum  is  given  by 
inf  72(t)  for  t  approaching  sup  t,  t  in  the  neighborhood. 
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If  g’(0)  =  0^  diverse  pathological  situations  can  occur.  It  is  expeditious  to  consider 
the  case  where  g(t)  =  ^{t),  (j){0)  f  0^  <j)  twice  continuously  differentiable  in  a  neighborhood 

of  t  =  0,  and  omit  discussion  of  any  pathological  situations.  (For  o'  <  1,  g*(0)  =  +oo,  ex¬ 
amination  shows  that  a  shock  appears  at  the  piston  at  t  =  0,  growing  continuously  from 
zero  strength.  )  The  interest  is  in  those  cases  where  a  >  1.  In  this  case. 


Ti(t)  -  T^(i)  = 


2 

7  +  1 


[^0 

1 

1  1 

[t^-l] 

a(^(t)  +  tf  (t)J 

r  5"^5)dt  - 


Hit)  _  I 


Since 


and 


C  =0(t),  =  o(t), 


lim 


1 


1 


a(^(t)  +  _r  ,^(0) 

there  exists  a  neighborhood  of  t  =  0  such  that 


Ti(t)  -  T2(t)  >  0; 

that  is,  in  this  neighborhood  is  to  the  left  of  Cq.  Both  T^(t)  and  T2(t)  tend 

to  +00  as  t  tends  to  zero.  Examination  of  shows  that  it  is  asymptotic  to  Cq. 

Three  final  observations  give  sufficient  information,  coupled  with  the  previous 
information,  to  find  the  infima  and  the  character  of  the  shock  which  will  occur  at  these 
various  points. 

If,  in  a  neighborhood,  g(t)  satisfies  equation  (17),  the  portion  of  the  envelope 
corresponding  to  this  neighborhood  is  a  point.  If  this  point  is  outside  the  region  of  influ¬ 
ence  of  any  other  shock  which  forms  in  the  flow,  an  outgoing  shock  appears  at  this  point 
with  jump  in  velocity  less  than  the  jump  of  g(t)  for  the  set  of  intersecting  characteristics. 
A  weaker  ingoing  shock  also  forms. 

If  g'(t)goes  to  zero  at  time  t^,  T2(t)  goes  to  +oo,  and  the  envlope  goes  to  infinity 
asymptotic  to  the  characteristic  corresponding  to  tj^,  x  =  x^(tp  +[cq  +  X.IJ:  g(tp](t  -  t^^). 

If,  while  g(t)  remains  continuous  at  t^,  g'(t)  is  discontinuous  at  t^  but  limits  of 
g'(t)  exist  to  the  right  and  left  of  t^,  the  value  of  T2(t)  is  discontinuous  at  tj(see  equation 
(9)),  as  is  that  of  x(t2).  Note,  however,  that  the  points  on  the  envelope  as  t  approaches 

from  the  left  and  right  are  both  on  the  characteristic  corresponding  to  t^. 

The  discontinuity  in  g'(t)  can  cause  the  envelope  to  shift  anywhere  along  the  corre¬ 
sponding  characteristic,  even  to  the  piston  front  (or  behind). 

The  spectrum  of  possibilities  arising  from  simple  waves  should  now  be  apparent. 

We  construct  now  a  g(t)  with  two  infima,  both  occurring  at  time  72*  Let  g(t)  be  linear. 
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g(0)  =  0,  up  to  the  time  t^.  Pick  g'(t^  +  0)  so  that  the  corresponding  point  on  the  en¬ 
velope  x/Tg)  ip  T^.  x(tj)  +[cq  +  V^g(ti)](T2  -  t^)  and  continue  g(t)  according  to 
equation  (16)  up  to  a  time  tg  <>2.  Set  g(t)  =  g(t2),  t  >  tg.  The  envelope  is  as  shown  in 
Fig.  7,  plus  the  point  at  infinity  along  the  characteristic  corresponding  to  tg. 


I 


At  time  Tg,  a  shock  forms  on  C^,  growing  from  zero  shock  strength.  Simultaneously 
at  Tg,  >^(t^)  +  Cq  +  T^t-^g(t^)(Tg  -  t^),  a  pair  of  shocks  with  positive  jump  in  velocity 
appears.  Note  that  a  contact  discontinuity  appears  with  the  latter  pair  of  shocks. 

Outside  the  range  of  influence  of  the  shocks,  the  simple  waves  remain  the  solution 
to  equations  (l). 


SUMMARY 

The  method  of  construction  which  is  presented  is  valid  for  all  simple  waves. 
Attention  is  focused  on  forward -facing  compression  waves.  Means  are  developed  for'  tire 
examination  of  the  infima  of  the  time  for  points  in  the  region  multiply  covered  by  the 
characteristics. 

The  author  has  made  cursory  examination  of  the  existence  theory  for  solutions  of 
equations  (1)  for  times  later  than  that  of  the  appearance  of  the  shock.  Adding  the  equa¬ 
tion 

X  =  XfTg)  +  ^  U  dt,  U  the  shock  velocity,  (19) 

■^2 

to  the  Volterra  integral  equations  for  continuous  flow  (requiring  the  Hugoniot  conditions 
hold  across  the  path,  velocity  and  pressure  continuous  behind  the  shock),  existence  of  a 
solution  with  a  shock  proceeds  readily  for  times  in  a  neighborhood  of  the  time  the  shock 
appears.  Moreover,  this  is  all  that  is  to  be  expected.  Return  to  the  simple  wave  whose 


envelope  appears  in  Fig.  7.  Had  g’(t,  +0)  been  such  that  its  point  on  the  envelope  corre¬ 
sponds  to  a  time  slightly  later  than  Tg,  the  corresponding  shocks  develop  at  this  slightly 
later  time. 

Only  by  adding  (for  a  forward-facing  compression  wave  with  "weak'*  piston  velocity) 
equation  (19)  with  U  defined  by 

“  lim  t)  -  p(X  -  e,  t)u(x  -  e,  t) 

U(X,t)-lun  p(X  +  e,t)-.p(X-c,t) 

(with  the  additional  requirements  above  on  density,  velocity,  .  ,  .  )  at  every  point  of  the 
fluid  does  it  seem  possible  to  discuss  existence  theory  in  any  global  sense  (all  time). 

Even  this  assumes  something  about  piston  strength.  That  is,  no  double  shocks,  as  in 
example  2,  occur.  In  this  later  case,  it  would  be  necessary  to  define  two  shock  velocities, 
U  and  U  ,  at  each  point,  and  accommodate  a  more  complex  algebraic  problem. 

Estimates  for  growth  of  shock  strength  in  terms  of  g(t)  are  possible,  and  will 
appear  in  a  subsequent  report. 


This  report  was  prepared  as  an  account  of  Government 
sponsored  work.  Neither  the  United  States,  nor  the  Com¬ 
mission,  nor  any  person  acting  on  behalf  of  the  Commission: 
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